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PACS 73.23.Hk - Coulomb blockade; single-electron tunneling 
PACS 71 . 10 .Pm - Fermions in reduced dimensions 

PACS 73 . 63 . -b - Electronic transport in nanoscale materials and structures 

Abstract. - We present a microscopic theory for interacting graphene armchair nanoribbon quan- 
tum dots. Long range interaction processes are responsible for Coulomb blockade and spin-charge 
separation. Short range ones, arising from the underlying honeycomb lattice of graphene smear 
the spin-charge separation and induce exchange correlations between bulk electrons - delocalized 
on the ribbon - and single electrons localized at the two ends. As a consequence, entangled end- 
bulk states where the bulk spin is no longer a conserved quantity occur. Entanglement 's signature 
is the occurrence of negative differential conductance effects in a fully symmetric set-up due to 
symmetry-forbidden transitions. 



The first successful separation of graphene [1], a sin- 
gle atomic layer of graphite, has resulted in intense theo- 
retical and experimental investigations on graphene-based 
structures [2], because of potential applications and fun- 
damental physics issues arising from the linear dispersion 
relation in the electronic band structure of graphene. 
In graphene nanostructures, confinement effects typical of 
mesoscopic systems and electron-electron interactions are 
expected to play a crucial role on the transport proper- 
l^ies. Indeed a tunable single-electron transistor has been 
demonstrated in a graphene island weakly coupled to leads 
['3]. Conductance quantization has been observed in 30nm 
v^ide ribbons [4], while an energy gap near the charge neu- 
trality point scaling with the inverse ribbon width was 
reported in [5]. Theoretical investigations [6,7] have at- 
tributed the existence of such a gap to Coulomb interac- 
tion effects. 

Confinement is also known to induce localized states at 
zig-zag boundaries [8], possessing a flat energy band and 
occuring in the mid of the gap. Those states have been 
analysed [9] under the assumption of a filled valence and 
an empty conduction band (half- filling) , taking into ac- 
count both Hubbard and long-ranged Coulomb interac- 
tion. There was a prediction of strong spin features in 
case of a low population of these midgap states. 

Above the half-filling regime, however, no detailed study 
on the interplay between longitudinal quantization ef- 
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Fig. 1: A graphene armchair nanoribbon single-electron tran- 
sistor. At the long sides, the lattice is terminated in armchair, 
at the small ends in zig-zag configuration. 



fects and Coulomb interactions in the spectrum of narrow 
nanoribbons exists at present. 

The purpose of this Letter is to derive a low energy the- 
ory of armchair nanoribbons (ACN) single-electron tran- 
sistors (SETs), see Fig. [U i.e., to investigate the conse- 
quences of confinement and interaction in narrow ACNs 
weakly coupled to leads. Short ACN have recently been 
synthesized [10]. We show that the long-range part of the 
Coulomb interaction is responsible for charging effects and 
spin-charge separation. Short-range processes, arising due 
to the presence of two atoms per unit cell in graphene as 
well as of localized end states, lead to exchange coupling. 
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Bulk-bulk short-range interactions have only a minor ef- 
fect on the energy spectrum. However, interactions be- 
tween end states localized at the narrow zig-zag ends of 
the stripe and bulk states smear the spin-charge separa- 
tion. Moreover, they cause an entanglement of end-bulk 
states with the same total spin. Hence, despite the weak 
spin-orbit coupling, the bulk spin is not a conserved quan- 
tity in ACNs. These states strongly influence the nonlin- 
ear transport. We predict the occurrence of negative dif- 
ferential conductance (NDC), due to symmetry-forbidden 
transitions between entangled states, in a fully symmetric 
setup. 

We proceed as follows: in the flrst part of this Letter we 
set up the interacting Hamiltonian of ACNs and derive 
their energy spectrum. In a second part transport in the 
single electron tunneling regime is investigated. 

Electron operator of a metallic ACN. — The car- 
bon atoms in graphene are arranged in a honeycomb lat- 
tice. There are two atoms per unit cell that deflne two 
different sublattices p = ±. Overlapping 2pz orbitals form 
valence and conduction 7r-bands that touch at the corner 
points of the flrst Brillouin zone, also called Dirac points, 
and determine the electronic properties at low energies. 
From now on we focus on the region of linear dispersion 
in the vicinity of the two inequivalent Dirac points, see 
Fig. Et, Kp = ^57=^6^, F = ±, where ao ^ 0.14 nm is 
the nearest neighbour distance. Then the 7r-electrons are 
described by Bloch waves 

=: ^ r?Fap(^) V^Fp(r,^), (1) 
p=± 

where Nl is the number of sites of the considered lattice, 
a = ± denotes the conduction/ valence band, and X^p(^) 
is the 2pz orbital on sublattice p at lattice site R, with 
f^R e M?. Furthermore = (hZx^i^y) is the wave vec- 
tor relative to the Dirac point Kp- Finally, the spinors 
VFa{i^) '= {VFa-{i^)^VFa-\-{i^)) fulflll the Dirac equation 
with a velocity vf = S.l ■ 10^ m/s. 

To describe ACNs boundary conditions have to be as- 
sumed. Following Ref. [11] we demand that the wavefunc- 
tion vanishes on sublattice p = — on the left end, y = 0, 
and on p = + on the right end, y = L. At the armchair 
edges the terminating atoms where the wave function is 
required to vanish are from both sublattices. The quanti- 
zation condition from the zigzag ends reads [11] 

that from the armchair edges is K-^ -\- = ^Ux, rix G Z. 
Eq. ([2]) supports the presence of extended states - real hZy 
- as well of localized states - purely imaginary K,y [8]. 

Let us flrst discuss the bulk states. Due to Ly Lx 
the longitudinal quantization condition yields subbands 




Fig. 2: a) Dispersion relation of a graphene stripe for real mo- 
menta. In the low energy regime, only subbands lying on the 
plane Kx = play a role due to the condition Lx Ly, with 
K the vector relative to the Dirac point (^-)- b) Confine- 
ment along the ACN length yields quantization of K,y. 

assigned to diflFerent Kx- From now on we focus on the low 
energy regime of metallic ACNs, where only the gapless 
subbands [k,x = 0, Fig. [2^)] are relevant. Eq. ([2]) yields 
then hZy = {fly + f^y ^ Z, Fig. [2)3). Bearing in mind 

Eq. (dl), we can finally express the states cp^y in terms of 
the sublattice wave functions (fFpj 

Fpr=± 

where r = ± denotes right /left moving waves. Up to a 
complex prefactor, the coefiftcients are /+r = ^, f-r = i- 

The quantization condition ([2]) also allows purely imag- 
inary hZy'. For each hZx = tIxTt/Lx > l/Ly, Ux G N 
there exist two imaginary solutions tvy{tvx)- Besides, due 
to Lx <C Ly, it holds to a very good approximation 
f^y{f^x) = ^ii^x- The corresponding ACN eigenstates can 
be chosen to live on one sublattice p = ± only: 

(^) = C{k.x) Y1 FifFpif, {Fkx, ipf^x))^ 

F 

where C{nx) is a normalization constant. The decay 
length of ^p^^ from one of the zigzag ends to the inte- 
rior is i^x^^ which is much shorter than the ribbon length. 
Hence end states are localized. From the graphene disper- 
sion relation it follows that the energy of the end states 
is zero. They will be unpopulated below half filling but 
as soon as the Dirac point is reached one electron will get 
trapped at each end. For small width ribbons the strongly 
localized character of the end states implies Coulomb ad- 
dition energies for a second electron on the same end by 
far exceeding the addition energy for the bulk states. Thus 
at low energies above the Dirac points both end states are 
populated with a single electron only. Introducing bulk 
and end electron annihilation operators c^j^iy^ d^jp,^^, the 
noninteracting Hamiltonian is 

Hq = hvF Kycl^^c^^y, (3) 
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because the end states have zero energy, and the fleld op- 
erator for an electron with spin a at position r is 



(4) 



= :'j/;^(r) 



The ID character of ACNs at low energies becomes evident 
by deflning the slowly varying electron operators ^pra{y) •= 



/2L, 



^'^ri^yV^^^^ g^(.]^ ^]^^^ obtain 



V^alO = jLy/2'^Ffpr(pFp{r)i^ra{y)^ 



where (pFp{r) :=(pFp{r^ i^ = 0). 



(5) 



Hamilton operator of the interacting ACN. — 

Including the relevant Coulomb interactions yields the 
total Hamiltonian 

He = Ho + Fe.b + V^b-b- (6) 
First, there is interaction between end and bulk states, 

K,x (Tcr'rr'p 

with y± = 0/L and with U^d{t — f') denoting the 3D 
Coulomb potential, the coupling constant 

Y,^^' II drdr'^l^{r)^F'+{r)UMr - r') 



e-b ^ / / 

FF' 



(7) 



For ACNs of width ranging from 5 to 25 nm, one flnds 
from numerical evaluation u^^_^ ^ '^e-b' with u^j^^Lxjs^ ~ 
0.55 nm, practically independent of Kx- 

Secondly, interaction between the extended bulk states, 

^b-b= E E ^fX^ 

is classifled by the scattering types Sr^Sa concerning 
band and spin, respectively, where one distinguishes be- 
tween forward (/^)-, back (6)-, and umklapp (u)- scat- 
tering. Denoting the scattering type by Sj we deflne 
[I]s,=f± ■■= [I,±I,±I,I], [I]b := [/,-/,/,-/] and [J]„ := 
[/, /, -/, -/], see also Fig. El With Eq. ® one finds 



b-b 



^{MsrA'^W} J I ^^^^ ^ 



X 4,Ay)4,,'{y')U[r]s^y,y')i^r,a'{y')i^r,.{y). (8) 

Hereby, the potential mediating the interactions is either 



^intra _ ^inter 







Fig. 3: The four different possibilities for scattering are forward 
(/^)-, back (6)-, and umklapp {u)- processes. As it can be seen 
from the sketch, those correspond to different relations between 
a certain index / of the states before and after the scattering 
event. 

where the ID potentials [/mtra/mter dggci-ibg interactions 
between electrons on the same/different sublattice [12]. 
While end-bulk scattering is completely short-ranged, the 
bulk-bulk interactions split into long-/short-ranged con- 
tributions {Sr=f^/Sr=u^b). The short-range bulk-bulk 
coupling constant is 



^b- 



(9) 



The long-ranged part of the interaction is diagonalizable 
by bosonization [13]. We flnd 

^0 + <ng = l^oN! + E(^<. + N^a) + H^os- (10) 

The flrst term of ([TQ|) . with Nc = J2a being the charge 
operator on the ACN, Wo = Wq=o with 

= ^ // ^y^y'^[r]f±iy^y')^^^iQy) ^^^iQy')^ 

accounts for Coulomb charging effects. The second term, 
where sq = hv/L is the level spacing, yields the ful- 
flllment of Pauli exclusion principle. Finally, i^bos — 
q>o ^jQ^lq^jQ ctccouuts for the bosonic excitations of 
the system, created/annihilated by the operators a^^ / 
ajq. The two channels j = c^s are associated to charge 
(c) and spin (s) excitations. The excitation energies are 

Esq = nqSo, Scq = UqEoy^l + 2Wq/eo with Tlq G N. 

Eigenstates of Hq + ^\^~^g I^L'^'^'^^)' where m 
characterizes the bosonic excitations, and the fermionic 
conflguration N = (N^^N^) deflnes the number of elec- 
trons in each spin band. Above half fllling exactly one 
electron occupies each end state and thus the end conflg- 
urations cr^,cr^ G {T^i}- 

These states can be used as basis to examine the ef- 
fect of V\^_]^ and Vq_]^ on the spectrum of an interacting 
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^ all short-range processes 
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Fig. 4: The spectrum of an ACN with A^c = 2n + 1 electrons. 
We chose a 7.8 nm width and 572 nm length, corresponding to 
a charging energy Wo = 2.3 meV, and to short-range bulk-bulk 
and end-bulk coupling constants u^_i, — 0.036 meV, Uq_^ — 
0.21 meV. End-bulk scattering i) mixes end and bulk states, 
ii) spoils the spin-charge separation. The inset zooms on the 
lowest lying 30 states. 



spin-modes are affected by short-ranged processes. The 
presence of end-states (a feature which is absent in car- 
bon armchair nanotubes [12]), however, smears out the 
energies within all plateaus (brown diamonds): It induces 
a mixing between excited states and groundstates of same 
total charge and spin, which widely lifts the degeneracy be- 
tween the various states. The inset of Fig. [H e.g., shows 
that among eight formerly degenerate groundstates, two 
get lowered and two get raised by a certain energy under 
the influence of the end-bulk interaction. We will come 
across this in more detail during the following analyis. 

Impact on transport. — In the remaining of this 
Letter we show how this entanglement is revealed in the 
peculiarities in the stability diagram of an ACN-SET. In 
the limit of weak coupling to the leads, we can assume 
that our total system, see also Fig. [H is described by the 
Hamiltonian 



ACN. For this purpose one needs to evaluate the corre- 
sponding matrix elements proportional to the short-range 
coupling constants ^^e-b' 0' %-b' ® • 
procedure follows similar lines as in [12] we refrain from 
reporting it here and discuss the main results. 
A diagonalization of the full Hamiltonian yields energy 
spectrum and eigenstates of the system including both 
long and short-range interactions. As those are spin pre- 
serving, it is clear that linear combinations must be formed 
of states with same spin-^^ component. Thereby, impor- 
tantly, the end spin degrees of freedom permit a mix- 
ture between states of different bulk spin conflgurations. 
This mechanism and its impacts will be illuminated in the 
course of the following sections. 

Spectrum of interacting ACNs. — Numerical cal- 
culation and diagonaliziton of the full ACN Hamiltonian 
including the 250 lowest lying states of a 572 nm x 7.8 nm 
ribbon leads to the spectrum found in Fig. [H For com- 
parison we also give the energies without the end-bulk 
interaction and for long-range interactions only. From Eq. 
([To]) it can be found that without short-range interactions 
(blue squares) , the energy cost for both a fermionic and a 
bosonic spin-like excitation amounts to sq. That is why in 
the spectrum discrete plateaus which are separated by this 
energy arise. The flrst charge-like bosonic mode can be ex- 
cited at an energy of about 2.l£o, which shows up in form 
of a small step towards the end of the third and all fol- 
lowing plateaus. Switching on the short-range bulk-bulk 
contributions (green disks) actually preserves this spin- 
charge separation: while the curve as a whole is shifted 
downwards in energy due to an exchange term (see inset 
of Fig. SI), all steps within the plateaus remain resolvable. 
In contrast to what is found for carbon nanotubes [12], 
there is only a very tiny additional lifting due to the bulk- 
bulk exchange, which cannot compare in magnitude with 
the spin-charge separation. The deeper reason is that, as it 
can be seen from an explicit calculation, only the bosonic 



H = Hq^ ^leads + - eQ^Vgate^^c 

with the ACN-Hamiltonian Hq given in Eq. (|6]). Fur- 
ther, i^ieads = T^igT^ai^q ' ^Maq^icrq. ^ith an- 
nihilating an electron in lead / of kinetic energy Cq and 
the chemical potential jii differs for the left and right 
contact by eV, with V the applied bias voltage. Next, 
Ht = Y^i^ Jd^r [Ti{r)^ljl{f)(j)ia{r) + h.c) descibes tunnel- 
ing between ACN and contacts, with tunneling coupling 
Ti{r) and V^cr(0 ACN bulk electron operator as given 
in Eq. g]), (t)icji{r) = Y.q^iqif)^icnq the lead electron op- 
erator with <j)iq{r) denoting the wave function of the con- 
tacts. Finally, the potential term describes the influence 
of a capacitively applied gate voltage (0 < a < 1). 
Due to the condition that the coupling between ACN and 
the contacts is weak, we can calculate the stationary cur- 
rent by solving a master equation for the reduced den- 
sity matrix to second order in the tunneling coupling. As 
this is a standard procedure, we refer to [14] for details 
about the method, and show in Fig. [5] numerical results 
for the differential conductance in the F-V^ate plane. In 
the numerical calculations an energy cutoff of 1.9eo above 
the groundstate was used, including any energetically al- 
lowed bosonic or fermionic excitation. One can clearly 
observe a two-fold electron periodicity, with small/large 
Coulomb diamonds corresponding to even, Nc = 2n, and 
odd, TVc = 2n + 1 electron fllling. A triplet of excitation 
lines is clearly visible in correspondence of the 2n 2n±l 
transition (Fig. O dashed red arrow). Moreover, NDC 
occurs as well, despite we considered a fully symmetric 
contact set-up (Fig. \S\ solid green arrow). To understand 
these features, it is necessary to consider the eigenstates of 
the fully interacting ACN in a minimal low-energy model. 

A minimal set of lowest lying states. — For the 

following we neglect short-range bulk-bulk processes as 
well as the bosonic excitations, as they do not qualitatively 
change the features we wish to describe. For even fllling, 
Nc = 2n, we consider those eigenstates |cr2,A^,cr^) : = 
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Fig. 5: Differential conductance of an ACN-SET. A triplet 
of states is split by the end-bulk interaction (dashed red ar- 
rows). The green solid arrows point towards regions where neg- 
ative differential conductance (NDC) is observable. We chose 
a temperature T= 116mK and tunneling coupling to the leads 
hrL = hrR=10~^ meV. All other parameters are as in Fig. 3. 

|cr^,7V,6,cr|) of Eq. ^ which have total spin 5 = 0, 
no bosonic excitations and up to one fermionic excitation. 
This means TV = (n, n) or TV = (n ± 1, n T 1), n G N. We 
introduce the notation (n, n) :=ti5 {n -\- l^n — 1) 
(n — l,n + l):=|| and get then four possible states, 



\a) :-- 



lTai,T), 



|5) :-- 
|c_) := 



U,Ti,T), 

U,TT,i>- 



The states \a)^\b) have the groundstate energy E^^'' 



£^2n 5 while the excited states \c±) have energy £^2n 
£^2n +^0- The mixing matrix elements, with Uq_]^ the end- 
bulk coupling constant, are Cl4_b)ac± = (^e-b)c±a = ^e-b' 
i^e-h)bc± = CKe-b)c±6 = "^e-b* Diagonalization yields: 



(/) 



Energy : 


Eigenstate {not normalized) 








^(|a)-|6)) + 


+ \c- 


-))=■■ 


|e2> 




\c+) - \c-) 










J^i\a)-\b)) + 


+ \c- 


-))=■■ 


|52) 




\a) + \b) 






l5i> 



2 V^2n 



aE. 



(0) 



iK-b)- 



where ^^a' - 2 V--2n ^ — 2n 

In total, the interaction has hardly lifted the degeneracies 
between the various states. However, symmetric and an- 
tisymmetric combinations of states |a), \b) and |c+), |c_) 
arise. The importance of this mixing becomes obvious 
when we look now at the states for the odd flllings. As 
we then necessarily have an unpaired spin, it is sufficient 
to consider merely the groundstates, i.e., TV = (n ± l,n) 
with energy £^2n+i total spin S = h/2. We introduce 
the notation, (n + 1, n) := t, — 1, := | and find the 
six states 



I«t) 



= IT,T,i), 
= li,T,T), 
= IT,i,T), 



I«t) 

\ci) 



= IT,i,i) 
= li,i,T), 
= li,T,i). 



The mixing matrix elements read G {T^ i})- 
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Fig. 6: Schematic explaining the mechanisms causing the NDC 
features (I) and (III) in Fig. \b\ Only states and transitions rel- 
evant for the NDCs are drawn. The crucial transition is marked 
by a big arrow head. (I) Opening of the channel |ti)^|e2) leads 
to a decay into the trapping state 1^2), depleting the transport 
channel \gi)^\ti). (Ill) Opening of the channel |t2)^|e2) de- 
pletes the transport channel 1^2) ^|t2). 



e: 



2n+l 



e: 



(0) 

2n+l 



Diagonalization yields: 

Eigenstate {not normalized) 
Its) := \a^) ^ \K) - V2c^ = 
M - \K) 
|a,) + \b,) + V2\c,) 



Vb- 

Energy 



p(0) 



\ts), 

\t2), 
\tl)- 



The excitation line triple. — Compared to the even 
fillings, the interaction induced lifting of the formerly de- 
generate 2n + 1 states is much more pronounced and seiz- 
able in the stability diagram of Fig. [5] in form of the triple 
of three parallel lines the dashed red arrow points to. The 
splitting has the expected value of V^Uq_]^. In detail, the 
lines mark transitions from the 2n groundstates \gi) and 
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Fig. 7: Occupation probability of the trapping state \g2) around 
the region exhibiting various NDC features. Their positions are 
marked according to Fig. [5] with labels (I)- (III). Notice that 
no numerically stable data can be obtained inside the Coulomb 
diamond. 

1^2) to the 2n + 1 states 1^2) and l^s). Hereby, the 
antisymmetric state 1^2), associated to the second line of 
the triple, is special, because it is the only one strongly 
connected to the 2n state |^2)- The first line of the triple 
is the 1^1 )^ 1^1 ) groundstate transition line. 

The NDC mechanism. - The NDC (I) highlighted 
by the solid green arrow marks the opening of the 2n+l ^ 
2n back-transition channel |ti)^|e2). The situation is 
sketched in Fig. [6l Once |e2) gets populated, from this 
excited 2n states the system can decay into any of the low- 
est lying 2n + 1 states, and in particular there is a chance 
to populate the antisymmetric state 1^2). This state is 
strongly connected to the 2n groundstate 1^2), which con- 
tains a large contribution of the antisymmetric combina- 
tion I a) — \h). But in the region where the NDC occurs, 
the forward channel 1^2) ^|^2) is not yet within the bias 
window such that \g2) serves as a trapping state. Fig. [71 
confirms this explanation: the population of the state \g2) 
is strongly enhanced in the concerned region where the 
back-transition |ti)^|e2) can take place, while the for- 
ward transition |^2)^|^2) is still forbidden. 
In a completely analog way, just involving instead of |e2) 
an excited 2n state with total spin h (not listed before), 
NDC (II) arises. 

The origin of NDC (III) is of different nature. It belongs 
to the back-transition 1^2)^ |e2), which is a weak chan- 
nel because \t2) is a purely antisymmetric state, while the 
antisymmetric contribution in |e2) is rather small. From 
time to time, nevertheless the transition will take place, 
and once it happens the system is unlikely to fall back to 
1^2), but will rather change to a symmetric 2n + 1 state. 
Thus the state \t2) is depleted, and with it the transport 
channel 1^2) ^1^2), which leads to NDC. The statement can 
also be verified from the plot of the occupation probabil- 
ity for 1^2), Fig. B a pronounced dark region of decreased 
population follows upon the NDC transition. 

Summary. - In conclusion, we focussed on small- 
width ACNs, and showed that the low energy properties 



are dominated by entangled bulk-end states. One major 
consequence is that the bulk spin is not conserved and 
that the symmetry of the entangled states generates trap- 
ping states and hence negative differential conductance. 

We acknowledge the support of the DFG under the pro- 
grams SFB 689 and GRK 638. 
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